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We investigate the complexation of a highly charged sphere with a long flexible polyelectrolyte, 
both negatively charged in salt free environment. Electroneutrality is insured by the presence of diva- 
lent counterions. Using molecular dynamics (MD) within the framework of the primitive model, we 
consider different Coulomb coupling regimes. At strong Coulomb coupling we find that the adsorbed 
chain is always confined to the colloidal surface but forms different conformations that depend on 
the linear charge density of the chain. A mechanism involving the polyelectrolyte overcharging is 
proposed to explain these structures. At intermediate Coulomb coupling, the chain conformation 
starts to become three-dimensional, and we observe multilayering of the highly charged chain while 
for lower charge density the chain wraps around the colloid. At weak Coulomb coupling, corre- 
sponding to an aqueous solvent, we still find like-charge complexation. In this latter case the chain 
conformation exhibits loops. 

PACS numbers: 61.20.Qg, 82.70.Dd, 87.10. +c 

I. INTRODUCTION 

The adsorpticaijaLpolyelectrolytes onto an oppositely charged spherical particle has txeen. gffPeri mmtal] y extensively 
studied recentlyum a-j-VaxkaJs authors have investigated this phenomenon theoreticallyH'Qu'BBEjOllda and by Monte 
Carlo simulationsolijIiSO'Ea. Nonetheless, much less is known concerning the complexation of a charged sphere with 
a like-charged polyelectrolyte. It is only .siery recently that we reported in a short communication such a phenomenon 
in the strong Coulomb coupling regimeO. From a theoretical point of view, the long range Coulomb interactions 
of these systems represents a formidable challenge, and especially the understanding of effective attraction of like 
charged bodies has attracted recent attention 

In this paper, we elaborate on the complexation between a sphere and a long flexible polyelectrolyte both negatively 
charged. While complexation in the strong coupling limit it is expected, we report new and rather unexpected 
chain conformations. We present MD simulation results without added salt but taking into account the counterions 
explicitly. Various Coulomb couplings as well as different polyelectrolyte linear charge densities are investigated. A 
detailed study of the ions (monomer and counterion) distribution is reported and mechanisms accounting for the 
different encountered complex structures are proposed. 



The paper is organized as follows. Section II contains details of our MD simulation model. Section III is devoted 



to the strong Coulomb coupling regime. Section IV is devoted to the intermediate Coulomb coupling regime. In Sec. 
0, we consider the like-charge complexation in the weak Coulomb coupling regime corresponding to water solvent. 

II. SIMULATION METHOD 

The iMD method employed here is based on the Langevin equation and is similar to that employed in previous 
studicso. Consider within the framework of the primitive model one spherical macroion characterized by a diameter 
d and a bare charge Qm = —Z^ie (where e is the elementary charge and Zm > 0) surrounded by an implicit solvent 
of relative dielectric permittivity e r . The polymer chain is made up of N m monomers of diameter er. Both ends of 
the chain are always charged and each 1/f monomer is charged so that the chain contains N cm — (N m — l)f + 1 
charged monomers. The monomer charge is q rn = —Z m e (with Z m > 0). The small counterions, assumed all identical, 
ensure global electroneutrality and have a diameter a and charge +Z c e (with Z c > 0). All these particles making up 
the system are confined in an impermeable spherical cell of radius R, and the spherical macroion is held fixed at the 
center of the cell. 

The equation of motion of any mobile particle (counterion or monomer) i reads 

(p'X- dv ■ 

m-^ = -ViU(n) ~ mi -± + W t -(i) , (1) 

where m is the mass particle (supposed identical for all mobile species), U is the total potential force and 7 is the 
friction coefficient. Friction and stochastic force are linked by the dissipation-fluctuation theorem < Wj(i)-Wj(i') >= 
&m'ykBT8ij5{t — t ). 
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Excluded volume interactions are introduced via a pure short range repulsive Lennard-Jones (LJ) potential given 

by 



U LJ (r) 




for r - r < r cut , , , 

for r - r > r cut , 



where tq — for the microion-microion interaction, ro = 7a for the macroion-microion interaction and r cut — 2 1/,6 cr 
is the cutoff radius. This leads to d = 2r$ + a = 15a, whereas the closest center-center distance of the microions to 
the spherical macroion is a = ro + a = 8a. The macroion volume fraction is defined as Jm = (a/R) 3 and is fixed to 
8 x 1CT 3 with R = 40er. 

The pair electrostatic interaction between any pair ij, where i and j denote either a macroion or a charged microion 
(counterion or charged monomer), reads 

Ucouiir) _ , ZiZj 



k B T 



Ib^, (3) 



where Ib — e 2 /AireoerksT is the Bjerrum length. Energies are measured in units of e = k B T with T = 298 K. 
Choosing a = 3.57 A rrequires that the Bjerrum length of water at room temperature equals 2a (7.14 A). In this 
work the macroion charge is fixed at Zm = 180. 

The polyelectrolyte chain connectivity is mO|deled by using a standard finitely extensible nonlinear elastic (FENE) 
potential in good solvent (see for example RefB3), which reads 



U F ENE{r) = 



— t^kRq In 
oo, 



for r < R Q , . . 

for r > Ro, K ' 



where k is the spring constant-like chosen as 1000fcsT/cr 2 and Ro — 1.5a. These values lead to an equilibrium bond 
length I = 0.8a. 

Typical simulation parameters are summarized in Table [j| The simulation runs are reported in Table H Each 
simulation run requires about 10 7 MD steps, and equilibrium is typically reached after 5 x 10 5 up to 3 x 10 6 steps. We 
normally performed between 5 x 10 6 and 10 7 MD steps to take measurements. We cover the whole range of strength 
of Coulomb electrostatic interaction from the strong coupling limit, which is more theoretical interest, to the weak 
coupling limit, which corresponds to an aqueous solvent. 



TABLE I: General data of the simulation model. Note that the temperature is used as energy scale for the simulations, 
parameters 

a = 3.57 A Lennard Jones length units 

T = 298K room temperature 

e = ksT Lennard Jones energy units 

Zm = 180 macroion valence 

Z m monomer valence 

Z c counterion valence 

7V C total number of counterions 

Ib Bjerrum length 
a — 8(7 macroion-counterion distance of closest approach 

R — 40a simulation cell radius 

fit = 8 x 1CP 3 macroion volume fraction 

k = lOOOfcT/cr 2 FENE spring constant 

Ro = 1.5cr FENE cutoff 

/ = 0.8a average bond length 

N m total number of monomers 

N cm number of charged monomers 

/ monomer charge fraction 

X PE — Z m ef/l polyelectrolyte linear charge density 
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TABLE II: Specification of the simulated systems. The chain radii of gyration R g bulk ' and R g comp ^ are given for an isolated 
chain (i.e., in the absence of the colloid) and for the complexation case (i.e., in the presence of the colloid) respectively. Lengths 
are in units of a. 
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III. STRONG COULOMB COUPLING 

First we look at the like-charge complexation in the strong Coulomb coupling regime. We choose the relative 
permittivity e r = 16, corresponding to Ib = lOc, and divalent microions (Z m = Z c = 2). Such a set of parameters is 
of special theoretical interest to study the influence of strong electrostatic correlations. 

A. Single charged object 

In this section we first study separately (i) the spherical macroion and (ii) the flexible polyelectrolyte in the presence 
of their surrounding neutralizing divalent counterions separately . This provides the reference states for the more 
complicated situation, where both of these two objects interact. 

1. Spherical macroion 

Consider an isolated macroion with its surrounding divalent counterions. A pertinent parameter to describe the 
Coulomb coupling for such highly asymmetric electrolyte solution (macroion and counterions) is the so-called "plasma" 
parameter T = Z^Ib/o, cc , where a cc (which will be determined below) is the [average distance (triangular lattice 
parameter in the ground state) between counterions lying on the macroion surfaceEH. For the strong Coulomb coupling 
considered we hnd r w 13 (with Zm — 180), and fpr-fipite. macroion volume fraction (here /m — 8 x 10 -3 ), all 
counterions lie in the vicinity of the macroion surfaceSE 3 iol_J. 

To characterize the counterion layer structure, we compute the counterion correlation function g(s) on the surface 
of the sphera__1l__l, defined as: 



c2 5W = (E 5 ( s '- s ^( s "-^))' ( 5 ) 

where c = Nc/Ana 2 is the surface counterion concentration (A^ = Zm/Z c being the number of counterions), and 
s = \s' — s"\ corresponds to the arc length on the sphere of radius a (center-distance of closest approach of macroion 
and counterion). Each counterion (located in the vicinity of the surface) is radially projected on the ("contact") 
shell around the macroion center of radius a — 8<r. Correlation functions are computed by averaging g(s) over 1000 
independent equilibrium configurations which are statistically uncorrelated. The pair distribution g(s) is normalized 
as follows 



c / 2TTsg(s)ds = (N c - 1). 
Jo 



(6) 
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FIG. 1: Counterion surface correlation function g(s) for the 90 divalent counterions (Zm = 180) in the case of an isolated 
macroion in strong Coulomb coupling (Ib = 10er). See Fig. H for a typical equilibrium snapshot. 




FIG. 2: Snapshot of the equilibrium counterion structure of an isolated macroion with a bare bare charge Zm = 180 [see Fig. 
[j]for the corresponding g(s)\. 

Because of the finite size and the topology of the sphere, g(s) has a cut-off at na (=25.1 a) and a zero value there. 
Therefore, at large values of s, o(s) cannot directly be compared to the correlation function in an infinite plane. 

Results are depicted in Fig. u\ for Zm = 180. The first peak 1 appears at about s = a cc w 3c whereas the second 
peak about s rj 6a and finally the third smalL-peak around s ~ 9a. This structure, which is highly correlated, is 
referred to as a strongly correlated liquid (SCL)E3I2j, but not yet a Wigner crystal. A typical equilibrium configuration 
is depicted in Fig. where one can see the local arrangement close to a triangular lattice. 



1 It is this value a cc fa 3cr which gives T as 13. 
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FIG. 3: Snapshot of the equilibrium conformation of an isolated polyelectrolyte chain made up of 256 monomers where all 
monomers are charged (/ = 1). Monomers are in white and counterions in dark grey. 



2. Polyelectrolyte chain 



Now we investigate an isolated polyelectrolyte chain together with its surrounding divalent counterions in bulk 
confined to the same spherical cell of radius R — 40a. We consider four monomer charge fractions (i. e. four linear 
charge densities) / = 1, 1/2, 1/3 and 1/5. The chain is made up of N m = 256 monomers (N m = 257 for / = 1/2). 
while Z m — Z c — 2, The polymer chain parameters are identical to those of the complexation case (see Table |[j]). The 
chain extension is characterized by its radius of gyration R g given by 



4 = ^(gc<-c«)'), m 

where tcm is the center of mass position of the chain. 

The corresponding values of R g can be found in Tabic [fl] (see R g bwlk ^ for runs A-D). The chain extension varies 
little with / and is roughly given by R g « 4er. For this strong Coulomb coupling, all counterions are "condensed" 
into the polyelectrolyte globule for all three linear charge densities i c o n s id ered. The strong counterion condensation 
induces a collapse of the chain, which is by now well understoodEZI'E3'r 9 lp H 31 l. A typical equilibrium chain conformation 
is shown in Fig. || for / = 1. As expected the structure is very compact and highly ordered. Very similar structures 
are obtained for / = 1/2, 1/3 and 1/5. 



B. Complexation 

We now investigate the complexation of a highly charged colloid with a long flexible polyelectrolyte, both nega- 
tively charged. Four different parameter combinations, denoted by run A, B, C and D, were investigated which are 
summarized in Table ||. Going from run A to D the polyelectrolyte charge fraction / decreases from 1 to 1/5. The 
contour length of the chain is much larger (N m l/d w 14 times) than the colloidal particle diameter. 



1. Observation of the complexation 



Figure |] shows typical equilibrium configurations of the colloid-polyelectrolyte complex. We notice that in all 
reported cases complexation occurs and the polyelectrolyte is completely adsorbed onto the colloidal surface, that is, 
in presence of a highly charged colloid, the polyelectrolyte conformation becomes quasi two-dimensional in contrast to 
the bulk case (compare Fig. [| with Fig. ||). However the structure of these resulting complexes depends strongly on 
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in red. 



4: Typical equilibrium configurations of the colloid-polyelectrolyte complex in strong Coulomb coupling (is = fOcr) for 
run ^4 (/ = 1), (b) run B (f — 1/2), (c) run C (/ = 1/3) and (d) run D (/ = 1/5). Monomers are in white and counterions 



the value of /. For the fully charged polyelectrolyte case [run A with / = 1 - see Fig. |(a)] the monomers are closely 
packed forming a two-dimensional compact Hamiltonian-walk with the condensed counterions on the polyelectrolyte. 
This structure consists of closed packed lines made from either counterions or monomers. When the linear change 
density is reduced [see Fig. |](b-d)], the complex structures are qualitatively different. In these cases the monomers 
are no longer closely packed. For run B [f = 1/2, Fig. ||(b)], the monomers spread more over the particle surface and 
the polymer partially wraps around the sphere exhibiting a quasi two-dimensional surface pearl-necklace structure. 
For run C [f = 1/3, Fig. |(c)] and run D [f = 1/5, Fig. f|(d)], the monomers spread entirely over the particle surface, 
and the chain wraps the colloidal particle leading to an almost isotropic distribution of the monomers around the 
spherical macroion. 

The like-charge complex formation is due to the strong counterion mediated correlations which are known to induce 
attractions in the strong Coulomb coupling regime. Basically, the charged species will try to order locally in a 
way which is compatible with the chain connectivity and the macroion surface constraints. We now quantify those 
observations and propose a simple mechanism to explain the observed conformations. 



2. Adsorption profile 

To quantify the adsorption of the monomers and counterions on the macroion particle surface, we analyze three 
quantities: (i) the ion radial distribution function ni(r), (ii) the ion fraction Pi(r) and (iii) the net fluid charge Qi(r) 
(omitting the macroion bare charge Zm), where r is the distance from the spherical macroion center. The radial ion 
distribution function ni(r) is normalized as follows 
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ni{r)Airr 2 dr = Ni with i~c,m; 



(8) 



where Ni is the total number of ions, and the subscripts c and m stand for counterion and monomer respectively. The 
reduced integrated ion radial density Pj(r) is linked to rij(r) via 



J r rii(r )A-Kr 2 dr 
Pi(r) = — — with i = c, m. 



(9) 



Note that for / < 1, neutral and charged monomers are all included in the quantities n m (r) and P m (r). The net fluid 
charge Q(r) is given by 



Q(r) 



Z c n c (r ) - Z m n m (r ) 



4:Trr 2 dr , 



(10) 



where we choose e = 1. 

Integrated distribution P(r), radial distribution n(r) and fluid net charge Q(r) profiles are depicted in Figs. ^|(a-c) 
respectively. Fig. ||(a) shows that all atoms for all runs are condensed within a distance of about 10<r from the colloid 
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FIG. 5: Ion adsorption profiles (runs A — D) as a function of the distance r from the macroion center, (a) Fraction P(r) and 
(b) radial density n(r) of counterions (thin lines) and monomers (thick lines) - (c) reduced net fluid charge Q{t)/Zm- 
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center and more than 80% of the monomers and counterions are within a distance of 9.3a from the colloid center 
corresponding roughly to two atomic layers [see also the radial distribution in Fig[|(b)]. Due to strong electrostatic 
attraction between the sphere and the counterions and strong electrostatic repulsion between the sphere and the 
charged monomers the first layer (r ~ a = 8a) is exclusively made up of counterions [see Figs ||(a-b)]. Note that 
the monomer depletion in this first counterion layer also concerns neutral monomers (runs B-D) and this effect is 
attributed to the chain connectivity. This means that the repulsion stemming from the charged monomers impose the 
polymer structure as long as / and the Coulomb coupling are sufficiently high, which is the case in the present study. 
The height of the first peak in the counterion n c (r) profile is almost independent on /. The second ion layer is mixed 
of monomers and counterions, however with a majority of monomers [see Figs ||(a-b)]. Indeed the first monomer peak 
in the n m (r) profile (located at r re 8.7a) and the second counterion peak in the n c (r) profile (located at r re 9.2a) are 
only separated by roughly Ar re 0.5a [see Fig |^(b)]. This leads to a medium position located at r re 9a corresponding 
to a bilayer thickness. The height of the second peak in the counterion n c (r) profile increases with increasing / [see 
Fig |(b)]. 

For / = 1 (run A), we observe a massive macroion charge inversion of more than 100% [i.e. Q(t)/Zm > 2] in the 
first layer as well as a strong charge oscillation [see Fig. [|(c)]. Upon reducing / (runs B — D) the macroion charge 
overcompensation decreases as well as the charge oscillation amplitude [see Fig. 0(c)] . This is due to the fact that 
upon reducing / less counterions are present and their correlations change. 

3. Poly electrolyte chain radius of gyration 

Next, we investigate the radius of gyration R g of the chain, i.e. Eq. (Q), in order to gain insight of the spreading of 
the monomers over the sphere. The results reported in Fig. || (for Ib = 10cr) show that R g increases with decreasing 
decreasing / which demonstrates that the spreading of the monomers over the macroion surface is enhanced by 
decreasing the polyelectrolyte charge density. The jump in R g is particularly large between / = 1 and / < 1. This is 
in agreement with the visual inspection of the chain conformations presented in Fig. ^. Moreover, the isotropic case 
(monomers fully spread over the particle) corresponding to R g re 8.7a ( 2 ) is already reached for / = 1/2 (run B). 
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FIG. 6: Radius of gyration of the polymer as a function of the polyelectrolyte charge fraction / for Ib ~ Wcr (runs A-D), 
l B = 4<7 (runs E-F) and l B = 2a (runs G-H). 
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2 This value corresponds to the the first monomer peak position in the n m (r) profile [see Figjd (b)] . 
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4- Surface counterion correlation function 



In this section we are interested in determining the structure of the "free" counterions which are not condensed 
onto the polyelectrolyte chain. Counterions are called "condensed" on the polyelectrolyte chain when they lie within 
a distance r c = 1.2a perpendicular to the chain (Fig. |^). All other counterions are called "free", although they are 
still adsorbed onto th e colloidal surface. To characterize the structure of the free counterions we proceed in the same 
way as in Sec. 



[II A 1. The surface free counterion correlation function gf ree (s) is now given by 



z2 free9free{s) 



]S(s'- Si )6(s"- Sj ) 



(11) 



where the sum in Eq. (|l] 
concentration, with Nf re 



) is restricted to the free counterions, and c/ ree = Nf ree /4ira 2 is the surface free counterion 
being the average number of free counterions. The normalization is obtained as follows 



Cfree I 2lTSg f ree (s)ds = {Nf ree - 1). 

/o 



(12) 



Results are depicted in Fig. g for runs A-D. The important result is that the first peak of gf re e(f) is located at 
the sa me posit ion (s « 3a) as in the "unperturbed" case of an isolated macroion (without polyelectrolyte) studied 
in Sec. |III A 1 (see Fig. 0). Although the second peak of 3/ ree (s) is less pronounced than in the "unperturbed" case 
(compare Fig. H with Fig. fil), the local order of the free counterion structure is still high as can be visually inspected 




FIG. 8: Free counterion (unbounded to the polyelectrolyte) surface correlation function gf r ee(s) for runs A-D. 
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on the snapshots sketched in Fig. [|. Thus the adsorbed chain only affects the counterion distribution significantly in 
its immediate neighborhood. 

5. Poly electrolyte overcharging 

We now show that the concept of polyelectrolyte overcharging can be used to explain the observed complex struc- 
tures. Let N c d be the number of counterions we consider as condensed onto the polyelectrolyte. Then the overcharging 
ratio xpe is defined as 



N r . 



XPE 



N r . 



(13) 



which is merely the ratio between the amount of the total condensed counterion charge and the polyelectrolyte bare 
charge. 

This "overcharging" can also be analytically predicted by the simple assumption that the presence of the polyelec- 
trolyte (with its counterions) does not affect the free counterion distribution 3 . Let us consider the bare charged chain 
plus its own neutralizing counterions as an uncharged object that gets overcharged by intercepting all counterions (of 
the macroion) whose center lie within a ribbon of width 2r c and area A r n, — 2r c N m l (Fig. 0). If c is the counterion 
(of the macroion) concentration then the theoretical overcharging ratio Xth is merely given by 



Xth 



1 



1 



N cm ZA-Ka? ' 



(14) 



and since the number of charged monomers N cm is given by N cm — (N m — 1)/ + 1, Eq. (nlh reduces for N m ^ 1 to 



with C 



2r c lZ M 



Xth ~ 1 + c/f, 



(15) 
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FIG. 9: Polyelectrolyte overcharge as a function of / (runs A-D). The dashed line corresponds to the theoretical prediction 
where Eq. (|l4[) was used. 



3 This mainly holds for the first correlational peak which is the most important for the present discussion, as we showed in the previous 
Sec. B 4 
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TABLE III: Polyelectrolyte overcharge xpe values as a function of / . 



Run 


1// 


Xpe - MD 


Xpe - Theory 


A 


1 


1.140 ±0.002 


1.21 


B 


2 


1.39 ±0.02 


1.43 


C 


3 


1.59 ±0.02 


1.64 


D 


5 


1.88 ±0.03 


2.06 



Results are presented in Fig. g and the corresponding values can be found in Table III. It indicates that in all 
cases overcharging occurs (i.e., xpe > 1), and that it increases with decreasing polyelectrolyte charge density. We 
have excellent agreement (less than 10% difference) between simulation results and our toy model [Eq. Jl4|)]. In turn 
it explains why xpe varies almost linearly with 1// in our simulations. 

The /-dependency of the complexation structure can be explained through the overcharging. For this we consider the 
overcharged polyelectrolyte as a dressed (or renormalized) chain [bare chain + counterionsLwith an (or renormalized) 
linear charge density X* PE = —(xpe — l)Ap_E that obviously has the opposite sign of AppE3. Similarly one can define 
the renormalized charge of a monomer as 

1m = -(XPE - l)<7m- (16) 

Using Eq. (|l6| ) and the results of Fig. ^ this shows that increases with increasing 1/ /. The overcharging leads to 
an effective local repulsion of the monomers, and subsequently to a bond stiffening of the chain. This_in turn explains 
why the chain expands with increasing 1// (see Fig. || and Fig. ^ for the corresponding structures )E3. 



IV. INTERMEDIATE COULOMB COUPLING 



In this section we are dealing with a higher dielectric constant (e r = 40), meaning that we consider weaker Coulomb 
coupling (lg — 4<t). Experimentally this could correspond to using alcohol as a solvent. We consider a set of three 
runs E — G with lg = 4<r (see Table [n]). Thus these systems are, up to a shorter Bjerrum length lg, identical to runs 
A — C. It will be helpful to start the discussion with the description of the observed complex microstructures. 



A. Complex microstructure 

Typical equilibrium macroion-polyelectrolyte complex structures are sketched in Fig. |l(| When the polyelectrolyte 
is fully charged (/ = 1, run E), Fig. [T^(a) shows again a strongly compact chain conformation. But in the present 
situation the chain does not spread on the macroion surface as it was the case in the strong Coulomb coupling [compare 
Fig. [Io|(a) with Fig. |J(a)]. In fact the conformation of the charged chain in presence of the macroion is very similar 
to the bulk conformation [compare [n](a) with Fig.J§] . 

By reducing the monomer charge fraction, Figs. [L0[(b-c) show that the complex microstructure for runs F, G is again 
qualitatively different from the fully charged case (run E). For these smaller polyelectrolyte linear charge densities 
(/ = l/2 and / = 1/3), the chain conformation is again almost wrapping around the colloid. For / = 1/2 (run F) the 
peal-necklace structure observed in the strong Coulomb coupling [see Fig. ||(b)] does not appear here, instead small 
loops appear [see Fig. |l0](b)]. For the smallest monomer charge fraction (/ = 1/3, run G) the monomers fully spread 
over the macroion surface and the conformation is not compact, and also small loops appear [see Fig. |l0|(c)]. Again 
it is observed that upon reducing the polymer charge density the chain expands, but this time it expands also into 
the radial direction away from the macroion. 

The forthcoming sections are devoted to study in more detail the monomer and counterion distributions. 



B. Adsorption profile 



Integrated distribution P(r), radial distribution n(r) and fluid net charge Q(r) profiles are given in Figs. fiT|(a-c) 
respectively. In a general manner, the fraction of monomers and counterions in the vicinity of the macroion surface 
is clearly smaller than the one obtained in the strong Coulomb coupling as expected [compare Fig. lfl(a) with Fig. 
|(a)]. Also the width in P(r) profile of adsorbed ions is enlarged with decreasing Coulomb couplingfcompare Fig. 
il|(a) with Fig. 0(a)]. These features show that the monomer adsorption is more diffuse (in the normal direction to 
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FIG. 10: Typical equilibrium configurations of the colloid-polyelectrolyte complex in moderate Coulomb coupling (Ib = 4cr) 
for (a) run E (/ = 1), (b) run F (f = 1/2) and (c) run G (/ = 1/3). Monomers are in white and counterions in red. 



the macroion sphere) as expected for a weaker Coulomb coupling. Concerning the ion density n(r) profile [see Fig. 
pjj(b)l) it is interesting to note that the height of the first peak in the monomer density n rn (r) profile is twice smaller 
for / = 1 than for / < 1, whereas the one from the counterion density profile is almost independent on /. In parallel, 
the macroion charge overcompensation as well as charge oscillation amplitudes are clearly reduced compared to the 
strong Coulomb regime [compare Fig. [ll](c) with Fig. |(c)] . 

For the fully charged polyelectrolyte case (/ = 1, run E) the n m (r) profile in Fig. [n](b) shows a strong second 
monomer peak and a weaker third one in agreement with the snapshot of Fig. [lO|(a) . The radial monomer ordering 
naturally goes along with a counterion ordering in anti-phase. In other words, multilayering of different chain segments 
occurs, but without strong adsorption of the macroion [compare Fig. |l0](a) with Fig. [|(a)]. Therefore we find three 
charge oscillations [see Fig. |ll](c)] against only two in the strong Coulomb coupling [see Fig. ||(c)]. 

For smaller linear charge density, Fig. [n](a) and Fig. |ll| (b) indicate that for runs F (/ = 1/2) and G (f = 1/3) 
the chain is almost fully adsorbed to the macroion surface without monomer chain multilayering [i. e., no appearance 
of monomer second peak in the n m (r) profile - see Fig. [Tl|(b)], However, the conformation is a little bit swollen 



probably due to the onset of loop formation, compare the snapshots in Fig. 10. As in the strong Coulomb coupling 
case, we find here only one charge oscillation for / = 1/2 and / = 1/3, and the amplitude of the reduced net fluid 
charge decreases with decreasing /. 
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C. Polyelectrolyte chain radius of gyration 

Results for the radius of gyration R g of the polymer chain are reported in Fig. || (for Ib = 4er). As for the high 
Coulomb coupling case with Ib — 10cr, R g increases with decreasing /. However, for / = 1, here we obtain R g ss 4.8tr 
which is clearly smaller than the value R g ~ 6.4<r obtained in run A (see Fig. |). This proves that for Ib — 4c the 
chain conformation stays as a globule, since we already found that for the same chain length and with Ib = 10cr the 
chain conformation was compact and two-dimensional. In the case under consideration (/ = 1, Ib = 4cr), the value 
of R g » 4.8(7 is almost identical to the one obtained in the bulk where R g — AA<r (see Table ||). 

Upon reducing / (run F and G), the chain is much more expanded and it is found that R g 9.4o\ Taking into 
account the fact the chain is still adsorbed for both systems (run F and G) as was found in the analysis of the 
adsorption profile (see Fig. pi] ), one deduces that the spreading of the monomers is very important as soon as / < 1. 
A comparison with the strong Coulomb coupling shows that values of R g for Ib = 4a are systematically larger than 
those of R g with Ib = 10cr (see Fig. || and Table ||), indicating again that the chain fluctuates more in the outward 
macroion radial direction at weaker Coulomb coupling. 
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TABLE IV: Polyelectrolyte overcharge Xpe values as a function of / (runs E — G). 



Run 


1// 


XPE 


E 


1 


1.043 ± 0.002 


F 


2 


1.089 ± 0.003 


G 


3 


1.065 ± 0.006 



D. Polyelectrolyte overcharging 



In order to check if the local polyelectrolyte overcharge is responsible for the expansion of the chain upon reducing / 
as was demonstrated in the strong Coulomb coupling for Ib = 10cr in Sec. Ill B 5| , we again consider the overc harging 
ratio xpe defined by Eq. ( |i~3| ) with the same condensation distance r c = 1.2c as was done for Ib = lOc in Sec. [II B I . 

Numerical values of Xpe are can be found in Table IV. It clearly shows that polyelectrolyte overcharge is negligible 
(i. e., xpe ~ 1) for the present Coulomb coupling regime whatever the value of /. Consequently one can not explain 
the expansion of the chain with increasing l/f with a polyelectrolyte overcharge mechanism. We will give clear and 
qualitative arguments in Sec. VI that account for these conformations. 



V. WEAK COULOMB COUPLING 



This part is devoted to aqueous solutions where the Bjerrum length is Ib = 2a = 7.14 A ccorresponding to the 
dielectric constant e r k, 80 of water. Such systems will be referred to as the weak Coulomb coupling regime. We have 
considered a set of three runs H — J (see Table [n]) identical to the previous ones but with a shorter Bjerrum length 
Ib = 2cr. Again it is helpful to start with the description of the observed complex microstructures. 



A. Complex microstructure 

Typical equilibrium macroion-polyelectrolyte complex structures can be found in Fig. |l2|. For all investigated cases 
(runs H — J), one finds that the polymer never adopts a "two-dimensional" conformation. A comparison of the bulk 
value of R g and the R g of the chain in the complexed situation (compare Table reveals that the chain conformation 
is only weakly affected by the macroion. For the fully charged polymer (/ = 1, run H), the conformation is again 
rather compact but without exhibiting a strong monomer-counterion ordering (within the polymeric aggregate) as it 
was the case for higher Coulomb coupling regimes [compare Fig. |l^(a) with Fig. [n](a) and Fig. 0(a)] . However we do 
have an effective macroion-polyelectrolyte attraction, and the dense monomer-counterion aggregate is adsorbed onto 
the colloidal surface. 

For / = 1/2 (run 7), the chain conformation is more expanded than for / = 1 [compare Fig. |l2|(b) with Fig. |i2|(a)]. 
Nevertheless we do have polymer adsorption with the formation of chain loops. Therefore even for couplings which 
are typical of aqueous systems, our simulations show that like-charge complexation can occur for / = 1 and / = 1/2 
with divalent ions (Z m = Z c = 2). 

For even lower linear charge density (/ = 1/3, run J), Fig. |l2|(c) shows polymer desorption from / = 1/2 to 
/ = 1/3. In the same time there is a certain "counterion release" for / = 1/3, meaning that not all counterions are 
in the vicinity of the highly charged objects (macroion and polyelectrolyte). We carefully checked that, with periodic 
boundary conditions, the same features qualitatively appear, namely chain desorption from / = 1/2 to / = 1/3. 



B. Adsorption profile 



Integrated distribution P(r), radial distribution n(r) and fluid net charge Q(r) profiles are depicted in Figs, (a-c) 
respectively. The ion fraction P(r) profiles show that for / = 1 (run H) and / = 1/2 (run I) almost all atoms lie 
within a distance a < r < 20a [i. e. P{r = 20a) « 1], corresponding to roughly one macroion diameter away from 
the colloidal surface [see Fig. 03(a)]. This is in contrast to what was previously found (with / = 1 and / = 1/2 ) 
at stronger Coulomb coupling regimes where almost all ions lie within a distance of a few monomer sizes from the 
macroion surface [compare Fig. P(a) with Fig. |(a) and Fig. 0(a)]. For / = 1/3 ( run J) only a very small fraction 
of monomers [P m (r = 10a) < 5%] lie in the vicinity of the macroion surface. In this latter situation, the counterion 
P c (r) profile indicates that a larger fraction of counterions float in the solution. Because the R g of the chain is 
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FIG. 12: Typical equilibrium configurations of the colloid-polyelectrolyte complex in weak Coulomb coupling (Ib = 2a) for 
(a) run H (/ — 1), (b) run I (f = 1/2) and (c) run J (f — 1/3). Monomers are in white and counterions in red. Snapshot (c) 
was obtained with periodic boundary conditions. 

very large and some chain monomers might be interacting with the cell boundary we performed for this situation a 
simulationjvhere we employed periodic boundary condition, and where the interactions were computed using the P3M 
algorithmEJ. With this run we found no chain monomers in the vicinity of the colloid surface, hence unambiguously 
found monomer desorption from the colloidal surface. 

Concerning the intra-chain monomer ordering for / = 1, Fig. [l^(b) shows that, although the chain conformation 
is relatively dense, there is only one main peak in the monomer n m (r) profile (the second peak is marginal). This 
proves that there is no strong intra-chain monomer ordering (in the normal direction to the macroion sphere) in this 
weak Coulomb coupling in contrast with our observations at Ib = 4<t [compare with Fig. p"i"|(b)]. Nevertheless a 
second counterion layer (third monomer layer) is build. The height h m of the monomer peak is identical (within the 
statistical uncertainty) for / = 1 and / = 1/2 and corresponds to h m f» 0.038/c 3 , whereas for / = 1/3 we have 
h m w 0.01/er 3 . 

As far as the net fluid charge Q(r) is concerned, Fig. [l3](c) shows that for f — I a weak charge oscillation appears 
with a marginal macroion overcharge compensation of 3%. For / = 1/2, the same marginal macroion overcharging 
occurs but without exhibiting charge oscillation. Finally, for / = 1/3 no overcharging appears and the the net charge 
increases monotonically. 
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FIG. 13: Ion adsorption profiles (runs H — J) as a function of the distance r from the macroion center, (a) Fraction P(r) and 
(b) radial density n(r) of counterions (thin lines) and monomers (thick lines) adsorbed on the spherical macroion. (c) Reduced 
net fluid charge Q(t)/Zm- 



C. Polyelectrolyte chain radius of gyration 

In Fig. H (for lg — 2a) the chain radius of gyration R g as function of / is plotted. It is found, that R g increases 
almost linearly with l/f. This result fits well with the scaling theory in this regime from Ref£3 where it is found 
that the chain extension shrinks proportionally with Ib, and we can assume for our purposes / oc Ib- Because the 
electrostatic interactions are weaker in the present case {Ib = 2cr), the ion pair (monomer - condensed counterion) 
attractions are weaker which results in a higher R g value (at fixed /) than in the strong Coulomb coupling regime. 



D. Monovalent case, / = 1 

Our last result concerns the monovalent case (Z m = Z c = 1 and / = 1) corresponding to run K (see Table ||. In 
this situation we have a strong macroion-polyelectrolyte repulsion, and one can not get like-charge complexation. 
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VI. CONCLUDING REMARKS 



We have carried out MD simulations to study the complexation of a charged colloid with a charged polyelectrolyte 
of the same charge for various Coulomb couplings lg and and varying charge fraction /. 

For lg = 10a we gave a reasoning for the observed conformation in terms of overcharging of the single chain. 
However this argument only worked for the largest coupling parameter. 

A complementary view of the observed conformations is to regard both macroions as being neutralized by their 
counterions The isolated chain would then collapse into a globule, and the colloid would be regularly covered by its 
counterions. By changing the Bjerrum length lg we change the correlations between the charges which lead in a first 
approximation to attractions of dipolar origin, and the attraction is roughly proportional to lg. On the other hand, 
by varying / we change the number of counterions of the chain, hence the number of available dipolcs. 

For / = 1 we can regard the colloidal particle as exerting only a perturbation of the chain complex (which has a 
higher density of dipoles) . For largest I g the strong attraction results in a flat disk, for weaker lg the disk swells back 
into the bulk structure, the globule. For smaller / the number of counterions of the polyelectrolyte and the colloid 
become comparable in number, and both macroions can equally compete for the counterions leading to a greater 
freedom of both macroions to "move" in their common counterion cloud. For the largest value of / g we have again 
the strongest dipolar attractions leading to a purely 2D conformation, where the chain is wrapped around the colloid. 
For smaller values of lg the conformations become more and more 3D-likc, but are still wrapping around the colloid. 
Of special interest is here the fact that even for a coupling strength which is typically for an aqueous solvent (lg = 7.1 
A) we find that like-charge complexation still occurs provided the linear charge density is sufficiently large. However 
the adsorption of the polymer chain onto the colloid is weaker than in the larger Coulomb coupling regimes. For 
/ = 1/2, we observe formation of large loops. 

For lg = 2a and / = 1/3 we find no complex, but end with two single macroions together with their counterion 
cloud which interact mainly with their bare charge, i.e. repulsive. 

Our parameters concerning aqueous solutions should be experimentally accessible, typically for small highly charged 
colloids (micelles and relatively short polyelectrolyte chains). 
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